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Abstract 



The Dirac equation for an electron in two spatial dimensions in the Coulomb 
and homogeneous magnetic fields is discussed. This is connected to the prob- 
lem of the two-dimensional hydrogen-like atom in the presence of external 



> 
OO 

o 

' magnetic field. For weak magnetic fields, the approximate energy values are 

o' 

C^) ■ obtained by semiclassical method. In the case with strong magnetic fields, 

we present the exact recursion relations that determine the coefficients of the 
series expansion of wave functions, the possible energies and the magnetic 
fields. It is found that analytic solutions are possible for a denumerably infi- 



nite set of magnetic field strengths. This system thus furnishes an example of 



^ ■ the so-called quasi-exactly solvable models. A distinctive feature in the Dirac 



case is that, depending on the strength of the Coulomb field, not all total an- 
gular momentum quantum number allow exact solutions with wavefunctions 
in reasonable polynomial forms. Solutions in the nonrelativistic limit with 
both attractive and repulsive Coulomb fields are briefly discussed by means 
of the method of factorization. 
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I. INTRODUCTION 



Planar nonrelativistic electron systems in a uniform magnetic field are fundamental quan- 
tum systems which have provided insights into many novel phenomena, such as the quan- 
tum Hall effect and the theory of anyons, particles obeying fractional statistics [0,0] • Planar 
electron systems with energy spectrum described by the Dirac Hamiltonian have also been 
studied as field-theoretical models for the quantum Hall effect and anyon theory ||. Related 
to these field-theoretical models are the recent interesting studies regarding the instability of 
the naive vacuum and spontaneous magnetization in (2+l)-dimensional quantum electrody- 
namics, which is induced by a bare Chern-Simons term Q. In view of these developments, 
it is essential to have a better understanding of the properties of planar Dirac particles in 
the presence of external electromagnetic fields. 

In || we studied exact solutions of planar Dirac equation in the presence of a strong 
Coulomb field, and the stability of the Dirac vacuum in a regulated Coulomb field. Quite 
recently, there appear interesting studies on the quantum spectrum of a two-dimensional 
hydrogen atom in a homogenous magnetic field 0,0. As is well known, hydrogen atom 
in a homogeneous magnetic field has attracted great interest in recent years because of its 
classical chaotic behavior and its rich quantum structures. The main result found in 0,^] is 
that, unlike the three-dimensional case, the two-dimensional Schrodinger equation || and 
the Klein-Gordon equation [f7| can be solved analytically for a denumerably infinite set of 
magnetic field strengths. The solutions cannot be expressed in terms of special functions 
(see also §). 

In this paper we discuss the motion of Dirac electron in two spatial dimensions in the 
Coulomb and homogeneous magnetic fields, and try to obtain exact solutions of a particular 
form. As in the case of the two-dimensional Schrodinger and the Klein-Gordon equation, by 
imposing a sufficient condition that guarantees normalizability of the wavefunctions (see the 
paragraph after eq. fl35|) ) , we can obtain the exact energy levels for a denumerably infinite set 
of magnetic fields. In the Dirac case, however, not all values of the total angular momentum 



2 



j allow exact solutions with the form of wavefunctions we assumed here. Solutions for the 
nonrelativistic limit of the Dirac equation in 2+1 dimensions are briefly discussed by means 
of the method of factorization. 

We emphasize that in this paper, by assuming an ansatz which guarantees normalizability 
of the wavefunction, only parts of the energy spectrum of the system are solved exactly. In 
particular, we do not obtain energy levels with magnitude below the mass value, which 
include the most interesting ground state solution. This is the same as in the Schrodinger 
and the Klein-Gordon case. All these three cases can therefore be considered as examples 
of the newly discovered quasi-exactly solvable models 0. In (3 + l)-dimension, no analytic 
solutions, even for parts of the spectrum, are possible so far. 

II. MOTION OF DIRAC ELECTRON IN THE COULOMB AND MAGNETIC 

FIELDS 

To describe an electron by the Dirac equation in 2+1 dimensions we need only three 
anticommuting 7 M — matrices. Hence, the Dirac algebra 

{V\ Y} = 2gT , gT = diag(l, -1, -1) (1) 

may be represented in terms of the Pauli matrices as 7 = cr 3 , 7 fc = ia^, or equivalently, 
the matrices (ai,a 2 ) = 7°(7 1 ,7 2 ) = (— cr 2 ,o"i) and (3 = 7 ||. Then the Dirac equation 
for an electron minimally coupled to an external electromagnetic field has the form (we set 
c = h = l) 

(idt-H D )V(t,r) = Q, (2) 

where 

H D = aP + (3m - eA° = o x P 2 - cr 2 P l + a 3 m - eA° (3) 

is the Dirac Hamiltonian, Pk = —id^ + eA^ is the operator of generalized momentum of the 
electron, A^ the vector potential of the external electromagnetic field, m the rest mass of 
the electron, and — e (e > 0) is its electric charge. The Dirac wave function 
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*(*,r) 



(4) 



/ 



is a two- component function (i.e. a 2-spinor). Here ipi(t,r) and ^(t, r) are the "large" and 
"small" components of the wave functions. 

We shall solve for both positive and negative energy solutions of the Dirac equation (0) 
and (Q) in an external Coulomb field and a constant homogeneous magnetic field B > 
along the z direction: 



A°(r) = Ze/r (e > 0), A x — -By/ 2, A y = Bx/2 . 



(5) 



We assume the wave functions to have the form 



*(t,x) 



exp(-iEt)i/ji(r,ip) 



(6) 



where E is the energy of the electron, and 



f(r)e 



1 1 ip 



(7) 



/ 



with integral number /. The function ipi(r,ip) is an eigenfunction of the conserved total 
angular momentum J z = L z + S z = —id /dip + cr 3 /2 with eigenvalue j — I + 1/2. One can 
of course consider wavefunctions which are eigenf unctions of J z with eigenvalues I — 1/2. 
These functions are of the forms of (^) with ipi given by 



V 



g(r)e 



(8) 



/ 



But ansatz (^|) is equivalent to ansatz ([7|) if one makes the change I — > I — 1. It should be 
reminded that Z is not a good quantum number. This is evident from the fact that the two 
components of ip\ depend on the integer I in an asymmetric way. Only the eigenvalues j of 
the conserved total angular momentum J z are physically meaningful. For definiteness, in 
the rest of this paper, all statements and conclusions, whenever angular momentum number 
I is mentioned, are made with reference to ansatz @ and (0). 



Substituting @ and (0) in (fj), and taking into account of the equations 



R ± = I — ± I- — 



5 



i 9 ei?r 



dr \r dip 2 



(9) 



we obtain 



df (I eBr\ „ / Za\ 

— -| I f 4- I R 4- T71 -I I 

dr \ r 2 



/+ (E + m+—j g = 



dg 
dr 



fl + l eBr\ /_ Za\ , n 

where a = e 2 = 1/137 is the fine structure constant. If we let 



(10) 



F{r) = \pr- f(r) , G(r) = \/r g(r) 



(11) 



eq. (|10"D becomes: 



dF _ f l + \ + e£r ' 
dr \ r 2 
rfG fl + h . eBr 
dr \ 



E + m + — )G = 



2 +^\G-{E-m+ —) F = . 



r 

Za 
r J 



(12) 
(13) 



By eliminating G in (12) and F in (|T3|), one can obtain the decoupled second order differential 
equations for F and G. At large distances, these equations have the asymptotic forms 
(neglecting r~ 2 terms): 



d 2 F 

dr 2 
d 2 G 



dr 2 



+ 
+ 



E 2 - m 2 - eB(l + 1) + - i(e£r) 2 

r 4 



F = 



„o o IEZol 1 , „ . 

E —m — eBl + e£r 2 

r 4 



G = 



(14) 
(15) 



The last term in these two equations, which is proportional to r 2 , may be viewed as the 
"effective confining potential" . 

The exact solutions and the energy eigenvalues with < E < m corresponding to 
stationary states of the Dirac equation (pi]) with B = were found in ||. The electron 
energy spectrum in the Coulomb field has the form 

-1/2 

E = in I , (16) 



K + yj(l + 1/2)2 _ {Za) 2y 



where the values of the quantum number n r are : n r = 0, 1, 2, . . ., if I > 0, and n r = 1, 2, 3, . . . 
if Z < 0. It is seen that 



E = myi - {2Za) 2 (17) 

for I = n r = 0, and i?o becomes zero at Za = 1/2, whereas in three spatial dimensions 
Eq equals zero at Za = 1. Thus, in two spatial dimensions the expression for the electron 
ground state energy in the Coulomb field of a point-charge Ze no longer has a physical 
meaning at Za = 1/2. It is worth noting that the corresponding solution of the Dirac 
equation oscillates near the point r — > 0. 

For weak magnetic field the wave functions and energy levels with E < m can be found 
from (|T2"|) and (|T3| ) in the semiclassical approximation. We look for solutions of this system 



in the standard form 

F(r) = A(r) exp(iS(r)) , G{r) = B(r) exp(tS(r)) . (18) 



Here A(r) and B(r) are slowly varying functions. Substituting ( |T8"D into (|T2"D and (|T3"D, we 
arrive at an ordinary differential equation for S(r) in the form 

2 / w \2 /I.I /r,\2 / „ D,vA 2 



T\ =Q = E 2 -m 2 - eB(l + 1/2) + *™L + ^ ~ " + ^ ^ 
CLT I r r 



The energy levels with E < m are defined by the formula 

Qdr = n \-J(l + 1/2) 2 - (Zaf + . = } , (20) 

1 X h w? + eB(l + 1/2) -E 2 \ * 



where r max and r min (r max > r min ) are roots of equation Q = 0. In obtaining (ppp, the term 
(eBr) 2 in Q has been dropped. If we require the energy spectrum to reduce to ( |16|) when 
£> = 0, we must equate the right-hand side of (|20D to 7rn r . As a result we obtain (for I ^ 0) 

1 -1/2 

(n r + y / (/ + 1/2)2 _(z«)2)2_ 
In the nonrelativistic approximation the energy spectrum takes the form 



eS / 1 

m H Z + - 

2m V 2 



(21) 



(ZoQ 2 m eg / 1\ 

non " 2(n r + |Z + l/2|) 2 + 2ml + 2> ' 
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Semiclassical motion of electron in the magnetic and Coulomb fields can be characterized 



by means of the so-called "magnetic length" l B = Jl/eB and the Bohr radius a-Q = 1/Zam 
of a hydrogen-like atom of charge Ze. When the magnetic field is weak so that Is S> o-b, 
or equivalently, B <C B cr = (Za) 2 m 2 /e, the energy spectrum is simply the spectrum of a 
hydrogen-like atom perturbed by a weak magnetic field. We obtain the Zeeman splitting 
of atomic spectrum depending linearly upon the magnetic field strength and the "magnetic 
quantum number" 1 + 1/2. 

In strong magnetic field the asymptotic solutions of F(r) and G(r) have the forms 
exp(— ar 2 /2) with a = eB/2 at large r, and r 7 with 



for small r. One must have Za < 1/2, otherwise the wave function will oscillate as r — > 
when I = and I = —1. In this paper we shall look for solutions of F(r) and G(r) which can 
be expressed as a product of the asymptotic solutions (for small and large r) and a series in 
the form 





(23) 



F(r) = r 7 exp(— ar 2 /2) V] 



(24) 



71=0 



G(r) = r 7 exp(— ar 2 /2) V] (3 n r 



(25) 



n=0 



with a 7^ , Po 7^ 0. Substituting (p^) and (p5[) into (0) and (|l"3"D , we obtain 




(26) 



(27) 



(71 + 7) 



(■ 




a n + Za/5„ + (E + m) (3 n -i - 2aa n - 2 



(n > 2) 



(28) 



from (|I2"D , and 




(29) 



(30) 



7 



from (0). 

Eq.fl26|) and ( p9| ) allow us to express /3 in terms of a in two forms: 

Za 



Za 



at 



a 



(31) 
(32) 



which are equivalent in view of the fact that 7 = J (I + 1/2) 2 — (Za) 2 . Solving (^) and 
(j30l) with n = 1 gives 

(l + l + l)(E-m) + (y + l + f)(E 



m) 



Oil 



Pi 



(2 7 + 1) (7 + I + \ 
2 (7 + Z) £ - m 



Za ao 



a . 



(33) 
(34) 



(27 + 1) 

From ( |3"DD one sees that /3 n (n > 1) are obtainable from a n and ct n _i. To determine the 
recursion relations for the a n , we simply eliminate (3 n and in ( p8|) by means of (j30l). 
This leads to (for n > 2): 



^ + 7 + I - (n 2 + 2n 7 ) a r 



+Za 



n + 7 + I - (E - m) + (n + 7 + Z + - ) (E + m 



+ U + 7 + Z + 



E 2 - m 2 - 2a [n + 7 + I 



1 



a n _ 2 = . 



(35) 



Following ||, we impose the sufficient condition that the series parts of F(r) and G(r) 
should terminate appropriately in order to guarantee normalizability of the eigenfunctions. 
It follows from (|35D that the solution of F(r) becomes a polynomial of degree [n — 1) if the 
series given by (|35|) terminates at a certain n when a n = a n+ i = 0, and a m = (m > n + 2) 
follow from d35|) . Then from ( |3~0|) we have /3 n+ i = /3 n +2 = . . . = 0. Thus in general the 
polynomial part of the function G(r) is of one degree higher than that of F. Now suppose 
we have calculated a n in terms of «o («o 7^ 0) from (|33"D and (^) in the form: 



a n = K(l, n, E, a, Z) a 



(36) 



Then two conditions that ensure a n = and a n+ \ = are 



K(l,n,E,a,Z) = (37) 

and 



E 2 -m 2 = 2a (n + j + l+ , n = l,2, ... (38) 

Since the right hand side of ( |3~8|) is always non-negative Q, we must have \E\ > m for the 
energy. We note here that, similar to the Schrodinger and the Klein-Gordon case, the 
adopted ansatz guarantees the normalizability of the wavefunction, but does not provide 
energy levels with magnitudes below \E\ = m. 

For any integer n, eqs. (|37|) and (|38"D give us a certain number of pairs (E,a) of energy 
E and the corresponding magnetic field B (or a) which would guarantee normalizability 
of the wave function. Thus only parts of the whole spectrum of the system are exactly 
solved. The system can therefore be considered as an example of the quasi-exactly solvable 
models defined in 0. In principle the possible values of E and a can be obtained by first 
expressing the a (or E) in (p7|) in terms of E (a) according to fl38|). This gives an algebraic 
equation in E (a) which can be solved for real E (a). The corresponding values of a (E) are 
then obtained from ([58|). In practice the task could be tedious. We shall consider only the 
simplest cases below, namely, those with n = 1, 2 and 3. In these cases, the solution of the 
pair (E, a) is unique for fixed Z and I. In general, for n > 3, there could exist several pairs of 
values (E, a) (cf. |||7j). Unlike the non-relativistic case, here negative energy solutions are 



possible. As in the case of the (3+l)-dimensional Dirac equation [|10j, the unfilled negative 
energy solutions are interpreted as positrons with positive energies. 

We mention once again that all the exact solutions presented below, including the restric- 
tions for the values of I (or more appropriately, the values of the conserved total quantum 



number j = I + 1/2), are obtained according to the ansatz (0), and (p^ ) and ( p5|) with poly- 
nomial parts. Exact solutions for the other parts of the energy spectrum, if at all possible, 
would require ansatz of different forms which are not known yet. 



1 For / > 0, this is obvious. For / < —1, one has —1/2 < 7 + / + 1 < 0, recalling that Za < 1/2. 
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1. n = 1. 



In this case we have «o 7^ and a n = (n > 1). From (^) one obtains the energies 

m , . 

E = - . 39 

2(7 + l + l) K J 

Eq.(pHl) with n = 1 then gives the corresponding values of magnetic fields a. This results 
show that, with the ansatz assumed here, solution with positive energy cannot be obtained 
with n — 1. Furthermore, the previously mentioned requirement that E < —m can only be 
met with / < 0. 



2. n 



We now consider the next case, in which cto, ol\ 7^ 0, and a n = (n > 2). This also 
implies (3 n 7^ (n = 0, 1,2) and (3 n = (n > 3). From (0), ( p5|) and (|33|) , we must solve 
the following set of coupled equations for the possible values of E and a: 

V 



E 2 - m 2 = 2a [ 2 + 7 + I + 



Za [(r + 1){E - m) + (T + 2)(£ + m)] ai + 2a (r + 2) a =0 
(2 7 + l)r«i + Za [T(E - m) + (V + 1){E + m)] a = . 



(40) 
(41) 
(42) 



Here r = 7 + / + l/2. From these equations one can check that E satisfies the quadratic 
equation 

2 7 + 1 . 



27 + 1 

(2r + i)(2r + 3)--^^r 



E + 4m(r + 1) E + m 2 



1 + 



(Za) 



. 



(43) 



This can be solved by the standard formula. One must be reminded of the constraint 
\E\ > m. For I > 0, we can obtain analytic solutions with both positive and negative 
energies. But when I < 0, analytic solutions can only be obtained for negative energy 
E < —m. Furthermore, it can be checked that \E\ is a monotonic decreasing (increasing) 
function of |/| (Za) at fixed Za (I). 



10 



For Za <C 1/2, i.e. for light hydrogen-like atoms, we can write down approximate 
expression for energy near the mass value, i.e. \E\ ~ m . We can obtain from ( (4"3| ) the 
approximate values of E: 

2(Za) 2 



E A 



ni 



1 + 



(2 7 + i)r 



(r + i)(r + 2) 



, l>0, 



(44) 



for positive energies, and 



-m 



2(Za) 



(T + l) 



Z > and I < , 



(45) 



(27 + 1) 

for negative energies (in fact, it can be checked from fl43|) that for I < 0, E is always close 
to — m for any Za < 1/2). 

When Za is close to Za = 1/2, we have \E\^> m for / > 0. In this case the energy E 
can be approximated by: 

. 2 (2r + l)(2r + 3) n_1/2 



E = ±m 



(Za) 



(46) 



(27 + i)r 

A consequence following from this formula is that, for each I > 0, there is a critical value of 
Z beyond which polynomial solution with n = 2 is impossible. The critical value of Z for 
each / is found by setting the expression in the square-root of ([46]) to zero. For / = and 
/ = 1, the critical values of Z are Za = 1/2.936 and 1/2.316, respectively. 

In the non-relativistic limit (see Sect. Ill), it is the upper, or the large, component /(r) 
of the Dirac wave function that reduces to the Schrodinger wave function. Hence, in order 
to compare with the results considered in 0, it would be appropriate to study the nodal 
structures of the function F(r) for positive energy solutions in the limit E ~ m. It is easy 
to see from ( |4"T| ) or (H2I) that in this limit, «o and «o have opposite signs. Thus F(r) has 
only one node in this limit, which is the same as in the Schrodinger case. 



3. n = 3. 



For the case of n = 3, exact solution of ([Tj]) and fl3"gp becomes much more tedious. Now 
the values of E and a are solved by the following coupled equations: 
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E 2 - m 2 = 2a (T + 3) , (47) 

Za [(r + 2)(E - m) + (r + 3)(E + m)] a 2 + 2a (V + 3) a x = , (48) 

4( 7 + l)(r + l)a 2 + Za [(T + 1)(E - m) + (T + 2){E + m)} a x +Aa (T + 2) a = , (49) 

(2 7 + l)r«i + Za [T(E - m) + (V + 1)(E + m)] a = . (50) 



In place of (|43"D we now have a cubic equation for the energy E. We shall not attempt 
to solve it here. It turns out that the equation satisfied by E can be reduced to quadratic 
ones without linear term in E in the low energy (E « m) and the high energy (E ^> m) 
limit, which correspond to small and large Z, respectively. The results are 



E 4 



m 



1 - 



2(Za) 2 (r + i)(r + 2)(r + 3) 
(2 7 + i)r(r + 2) + 2( 7 + i)(r + 1) 2 



-1/2 



(51) 



and 



for \E\ ~ m, and 



-m 



E = ±m 



1 - 



2(Za) 2 (r + i)(r + 2)(r + 3) 
(2 7 + i)(r + 2) 2 + 2(7 + i)(r + i)(r + 3) 



2 {Zaf (r + i/2)(r + 3/2)(r + 5/2)(r + 3) 



1-1/2 



(52) 



-1/2 



(53) 



(2 7 + i)r(r + 5/2)(r + 2) + 2(7 + i)(r + i/2)(r + i)(r + 3)_ 

for \E\ ^> m. The corresponding values of the magnetic field are obtained by substituting 



(|1|), (0), or © into For I = -1, eq.fl53|) is real only for 1/2.65 < Za < 1/2. 

As in the n = 2 case, we shall also investigate the nodal structures of the function F(r) 
for positive energy solutions in the limit E ~ m. The zeros of the polynomial part of F(r) 
is given by 



1 

2 



41^ 



(54) 



,a 2 y V \«2/ V«2- 
Note that physical solutions of r , if exist, must be non-negative. In the limit E ~ m, 
eq.(f48"D and (0) give approximately 



a 'i 

«2 
«2 



a 

(2 7 + i)r «i 

2EZ«(r + 1) a 2 
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(55) 
(56) 



We see from (|55|) that ai/«2 < in this limit. 

For negative Z < 0, which implies — 1/2<T<0, we also have ao/«2 < 0. Eq.fl5l|) then 
implies that there is only one positive zero of F(r). Hence the wave function has only one 
node for / < 0. 

When Z > 0, we have T > 0, and hence «o/ a 2 > 0. It can be checked from (pETp, (pT|). 

that (ai/a^) 2 > 4(« /«2)- Thus F(r) has two positive zeros. This is also 
consistent with the results presented in || for the Schrodinger case (see also the last part of 
the following Section). 



and (56 



III. NON-RELATIVISTIC LIMIT AND METHOD OF FACTORIZATION 



The electron in 2+1 dimensions in the nonrelativistic approximation is described by 
one-component wave function. This can easily be shown in full analogy with the (3+1)- 
dimensional case. Let us represent \l/ in the form 



^ = exp(— imt) 



(57) 



and substitute fl5"T|) into (0). This results in, to the first order in 1/c, the following 
Schrodinger-type equation (instead of the Schrodinger-Pauli equation in 3+1 dimensions): 

.dip 



dt 



2m 2m r I 



(58) 



where, as before, Pu 



-id^ + eA^ denote the generalized momentum operators. The 



term eB/2m in (p8|) indicates that the electron has gyromagnetic factor g = 2 as in the 
(3 + l)-dimensional case [fLOR . 

One can now proceed in the same manner as in the Dirac case to solve for the possible 
energies and magnetic fields. We shall not repeat it here. More simply, we make use of the 
fact that eq.([58"D differs from the Schrodinger equation discussed in || only by the positive 
spin correction term lol = eB/2m, which is the Larmor frequency. We thus conclude that 
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the denumerably infinite set of magnetic field strengths obtained in || are still intact, but 
the corresponding values of the possible energies are all shifted by an amount Ul, i.e. 



E = co L (n + l + l + 



(59) 



Simply put, the quantum number n in || is changed to n + 1. 

Let us note here that the energies and magnetic fields in this case may also be found by 
means of a method closely resembling the method of factorization in nonrelativistic quantum 
mechanics. We shall discuss this method briefly below. Both the attractive and repulsive 
Coulomb interactions will be considered, since planar two electron systems in strong external 
homogeneous magnetic field (perpendicular to the plane in which the electrons is located) 
are also of considerable interest for the understanding of the fractional quantum Hall effect. 
Let us assume 



V(t,x) 



1 



^7T 



exp(—iEt + il(p)r^ exp(—ar 2 /2)Q(r 



(60) 



where Q is a polynomial, and a = eB/2 as defined before. Substituting ( |5"UD into flSlf), we 
have 



* 2 + (*l_ x \4. + U!l 



dx 



x 



dx 



x 



Q{x) = 



(61) 



Here x = r/ls, Ib = l/ve5, 7 = \l\ + 1/2, b = 2m\Z\alB = \Z\a^j2m/ui L} and e = 
E/lol — (2 + I + The upper (lower) sign in (|61~1) corresponds to the case of attractive 
(repulsive) Coulomb interaction. This will be assumed throughout the rest of the paper. 

It is seen that the problem of finding spectrum for ([H]) is equivalent to determining the 
eigenvalues of the operator 



H 



dx 2 



2 7 \ d b 

— x hr T_ • 

x J dx x 



(62) 



We want to factorize the operator (|62[) in the form 



H = a + a + p, 



(63) 
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where the quantum numbers p are related to the eigenvalues of (|6lD by p = e. The eigen- 
functions of the operator H at p = must satisfy the equation 

aip = . (64) 

s 

Suppose polynomial solutions exist for fl6jj) , say Q = II ~ x k), where Xk are the zeros of 

fc=i 

Q, and s is the degree of Q. Then the operator a must have the form 

8 s 1 

a = #■ - £ - (65) 

era ^ a; - x fc 

and the operator a + has the form 

« + = -#-- + *-E^— • (66) 



Substituting and (HBj) into ( |5"3"D and then comparing the result with Q5"^), we obtain 
the following set of equations for the zeros Xk (the so-called Bethe ansatz equations II): 

^-x fc -2^^^ = 0, k = l,...,s , (67) 

as well as the two relations: 

±6 = 2 T £X 1 , s = p. (68) 

it=i 

Summing all the s equations in ( |67| ) enables us to rewrite the first relation in ([68]) as 

±&=f>*. (69) 
fc=i 

From these formulas we can find the simplest solutions as well as the values of energy and 
magnetic field strength. The second relation in ( |68|) gives E = ul(2 + s + 1 + \l\), which is 
the same as in (|59"D noting that n = s + 1. 

For s = 1, 2 the zeros Xk and the values of the parameter b for which solutions in terms 
of polynomial of the corresponding degrees exist can easily be found from (|67|) and (|69"D in 
the form 



8 = 1, x 1 = ±J2\l\ + l, b = J2\l 



8 = 2, x x = (2|Z| + l)/x 2 , x 2 = ±(1 + y/4\l\+3)/y/2 , b = ^2(4|/| + 3) . (70) 
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From (|70|) and the definition of b one has the corresponding values of magnetic field strengths 

(Za) 2 

— & -2, (71) 



as well as the energies 



The corresponding polynomials are 

Qi = x — xi=x=f&, 

2 

^2= n(x-x fc ) =x 2 T^ + 2|/| + l . (73) 

k=l 

The wave functions are described by flSOf). For s — 1,2 for the repulsive Coulomb field the 
wave functions do not have nodes (for \l\ = 0,1), i.e. the states described by them are 
ground states, while for the attractive Coulomb field the wave function for s = 1 has one 
node (first excited state) and the wave function for s = 2 has two nodes (second excited 
state). 



IV. CONCLUSIONS 

In this paper we consider solutions of the Dirac equation in two spatial dimensions in the 
Coulomb and homogeneous magnetic fields. It is shown by using semiclassical approximation 
that for weak magnetic fields all discrete energy eigenvalues are negative levels of a hydrogen- 
like atom perturbed by the magnetic field. For large magnetic fields, analytic solutions of 
the Dirac equation are possible for a denumerably infinite set of magnetic field strengths, 
if the two components of the wave function are assumed to have the forms ( PH) and ( |25| ) 
with terminating polynomial parts. Such forms will guarantee normalizability of the wave 
functions. We present the exact recursion relations that determine the coefficients of the 
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series expansion for solutions of the Dirac equation, the possible energies and the magnetic 
fields. Exact and/or approximate expressions of the energy are explicitly given for the three 
simplest cases. For low positive energy solutions, we also investigate the nodal structures 
of the large components of the Dirac wave functions, and find that they are the same as 
in the Schrodinger case. We emphasize that, by assuming a sufficient condition on the 
wavefunction that guarantees normalizability, only parts of the energy spectrum of this 
system are exactly solved for. In this sense the system can be considered a quasi-exactly 
solvable model as defined in |§. As in the Schrodinger and the Klein-Gordon case, energy 
levels with magnitude below the mass value, which include the most interesting ground state 
solution, cannot be obtained by our ansatz. For the corresponding case in (3 + l)-dimension, 
no analytic solutions, even for parts of the spectrum, are possible. 
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